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Irreversibility is one of the most intriguing concepts in physics. While microscopic physical laws are perfectly 
reversible, macroscopic average behavior has a preferred direction of time. According to the second law of 
thermodynamics, this arrow of time is associated with a positive mean entropy production. Using a nuclear 
magnetic resonance setup, we measure the nonequilibrium entropy produced in an isolated spin-1/2 system 
following fast quenches of an external magnetic field and experimentally demonstrate that it is equal to the 
entropic distance, expressed by the Kullback-Leibler divergence, between a microscopic process and its time- 
reverse. Our result addresses the concept of irreversibility from a microscopic quantum standpoint. 


The microscopic laws of classical and quantum mechanics 
are time symmetric, and hence reversible. However, paradox¬ 
ically, macroscopic phenomena are not time-reversal invari¬ 
ant mm. This fundamental asymmetry defines a preferred 
direction of time that is characterized by a mean entropy pro¬ 
duction. Regardless of the details and nature of the evolution 
at hand, such entropy production is bound to be positive by 
the second law of thermodynamics 0. Since its introduction 
by Eddington in 1927 |4|, the thermodynamic arrow of time 
has not been tested experimentally at the level of a quantum 
system. 

Introduced by Clausius in the form of an uncompensated 
heat, the importance of the entropy production in nonequilib¬ 
rium statistical physics has been recognized by Onsager and 
further developed by Meixner, de Groot and Prigogine (5). 
Defined as E = /3(W — A F), for a system at constant inverse 
temperature (3 = 1 /(fc^T), where W is the total work done 
on the system and A F the free energy difference (kp is the 
Boltzmann constant), it plays an essential role in the evalu¬ 
ation of the efficiency of thermal machines, from molecular 
motors to car engines 0 . 

Starting with Boltzmann’s work on the so-called 77- 
theorem, the quest for a general microscopic expression for 
the entropy production, especially far from equilibrium, has 
been a challenge for more than a century 0. In the last years, 
formulas for the entropy production and entropy production 
rate in terms of the microscopic density operator p of a sys¬ 
tem have been obtained for relaxation m, transport ta, and 
driven processes in closed and open quantum systems EH3. 
At the same time, the recent development of fluctuation the¬ 
orems iflQl flTTl has led to a sharpening of the formulation of 
the second law. Regardless of the size of a system, the arrow 
of time originates from the combination of an explicit time- 
dependence of the Hamiltonian of the system and the specific 
choice of an initial equilibrium state. While the first ingredi¬ 
ent breaks time homogeneity (thus inducing the emergence of 
an arrow of time), the second specifies its direction lfl2l . 

In small systems, thermal and quantum fluctuations are 


both significant, and fluctuation theorems quantify the occur¬ 
rence of negative entropy production events during individual 
processes CG). In particular, the average entropy production 
(E) for evolution in a time window r has been related to the 
Kullback-Leibler relative entropy between states p v t and p^_ t 
along the forward and backward (i.e. time reversed) dynamics 
lfl4l-16l (see Fig. |T]). Explicitly 

(S) = S (p F t || p*_ t ) = tr [p F t (In p F t - In p F _ t )] . (1) 

The above equation quantifies irreversibility at the micro¬ 
scopic quantum level and for the most general dynamical pro¬ 
cess responsible for the evolution of a driven closed system. A 
process is thus reversible, (E) = 0, if forward and backward 
microscopic dynamics are indistinguishable. Nonequilibrium 
entropy production and its fluctuations have been measured 
in various classical systems, ranging from biomolecules Ifl7l 


Forward 



Figure 1. A quantum system (with Hamiltonian Ho) is initially 
prepared in a thermal state p e f at inverse temperature f3. It is driven 
by a fast quench into the nonequilibrium state p ¥ r along a forward 
protocol described by the Hamiltonian 77?. In the backward process, 
the system starts in the equilibrium state pf corresponding to the 
final Hamiltonian 77? and is driven by the time-reversed Hamiltonian 
77? = —H ¥ - t to p?. The entropy production (E) at time t is given 
by the Kullback-Leibler divergence between forward and backward 
states p ¥ t and p?_ t [cf. Eq. ([l}]. 
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and colloidal particles da to levitated nanoparticles ED (see 
Refs. l20ll2T1 for a review). Evidence of time asymmetry has 
been further observed in a driven classical Brownian particle 
and its electrical counterpart [ 221. However, quantum experi¬ 
ments have remained elusive so far, owing to the difficulty to 
measure thermodynamic quantities in the quantum regime. To 
date, Eq. G has thus never been tested. 

In this Letter, we use a nuclear magnetic resonance (NMR) 
setup to provide a clear-cut assessment of Eq. G where (E) 
and S (pt\\ Pr-t) are tested and evaluated independently. Our 
methodological approach is founded on the reconstruction of 
the statistics of work and entropy, following a non equilibrium 
process implemented on a two-level system, therefore assess¬ 
ing the emergence of irreversibility (and the associated arrow 
of time) starting from a genuine microscopic scale. 

We consider a liquid-state sample of chloroform, and en¬ 
code our system in the qubit embodied by the nuclear spin of 
13 C [23-25]. The sample is placed in the presence of a lon¬ 
gitudinal static magnetic field (whose direction is taken to be 
along the positive z axis) with strong intensity, Bq ~ 11.75 T. 
The 13 C nuclear spin precesses around Bq with Larmor fre¬ 
quency ccc/27t « 125 MHz. We control the system magneti¬ 
zation through rf-field pulses in the transverse (x and y) direc¬ 
tion EJ. The initial thermal state pg q of the 13 C nuclear spin 
(at inverse temperature /3) is prepared by suitable sequences 
of transversal radio-frequency (rf)-field and longitudinal field- 
gradient pulses. We use the value ksT/h = 1.56 ± 0.07 kHz 
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Figure 2. We show a section of the magnetometer employed in 
our NMR experiment. A superconducting magnet, which produces 
a high intensity magnetic field (Bo) in the longitudinal direction, is 
immersed in liquid He, surrounded by liquid N in another vacuum 
separated chamber, in a thermally shielded vessel. The liquid sample 
(inside a 5mm glass tube) is placed at the center of the magnet within 
the rf-coil of the probe head. A digital electronic time-modulated- 
pulse induces a transverse rf-field (B±) that drives the 13 C nuclear 
spins out of equilibrium. In the forward (backward) protocol, the 
rf-coil can perform (can retrieve) work on (produced by) the nuclear 
spin sample. Depending on the speed of the modulation of the driv¬ 
ing field, irreversible entropy is produced. The sketch is not in scale 
and has been stripped of unnecessary technical details of the appara¬ 
tus. 
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Figure 3. (a) [(b)] Evolution of the Bloch vector of the forward 

[backward] spin-1/2 state p ¥ t [p®_ t ] during a quench of the transverse 
magnetic field, obtained via quantum state tomography. A sampling 
of 21 intermediate steps has been used. The initial magnetization 
(gray arrow) is parallel to the external rf-field, aligned along positive 
x [y] axis for the forward [backward] process. The final state is rep¬ 
resented as a red [blue] arrow, (c) Polar projection (indicating only 
the magnetization direction) of the Bloch sphere with the trajectories 
of the spin. Green lines represent the path followed in a quasistatic 
(r —>> oo) process. 


(corresponding to an effective temperature of T ~ 75 =b 3nK) 
throughout the experiment for the initial 13 C nuclear spin ther¬ 
mal states. A sketch of the experimental setup is provided in 
Fig-S 

The system is driven out of equilibrium to the state p ¥ r by 
a fast quench of its Hamiltonian (denoted as H v t in this for¬ 
ward process) lasting a time r. We experimentally realize this 
quench by a transverse time-modulated rf field set at the fre¬ 
quency of the nuclear spin. In a rotating frame at the spin Lar¬ 
mor frequency, the Hamiltonian regulating the forward pro¬ 
cess is 


= 27rhiy (t) [cr£ cos + dy sin (2) 

with cr^ yz the Pauli spin operators, f(t) = nt/(2r ), and 
v(t) = vo (1 — t/r ) + v T t/r the (linear) modulation of the 
rf-field frequency over time r, from value vo — 1-0 kHz to 
v T = 1.8 kHz. With these definitions, the initial thermal state 
of the 13 C system is p^ = exp (—/3'Hq) /Zo, where Zo is the 
partition function at time t = 0. 
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Figure 4. (a) Black dots represent the measured negative and positive values of the entropy production E of the spin-1/2 system after a quench 
of the transverse magnetic field of duration r = 100 fis. The mean entropy production (red dashed line) is positive, in agreement with the 
second law. (b) Average entropy production (E) (dashed lines) evaluated through the probability distribution P(E), and Kullback-Leibler 
divergence S (pt 11 Pr-t) between forward and backward states p ¥ t and p B _ t (dots), reconstructed through the tomographic measurements, 
as a function of time for three different quench durations r = 100 ps (blue), 500 ps (green), and 700 fi s (red). Good agreement (within 
experimental uncertainties represented respectively by error bars and shadowed regions) between the two quantities is observed, quantifying 
the arrow of time Eq. (1). (c) Mean entropy production (E) (blue dots) and its linear response approximation (Elr) (green dots) as a function 
of the quench time. The difference between the two values, especially for fast quenches (small r), suggests the systematic deviation of the 
experiment from the linear response regime. The dashed lines represent the results of numerical simulations. 


In order to reconstruct the work and entropy production 
statistics of the 13 C quenched dynamics, we use the method 
proposed in Refs. l26l - [28l and illustrated in detail in the Ap¬ 
pendix. We consider an ancillary qubit embodied by the 
nuclear spin of our sample (Larmor frequency cch/27t « 
500 MHz) and exploit the natural scalar coupling between 
the X H and 13 C nuclear spins to implement the interferometer 
needed to reconstruct the statistics of the work done by 13 C 
following the quench (cf. Appendix). The method assumes a 
unitary dynamics of both system and ancilla, a condition that 
is met with excellent accuracy in our experiment. In fact, the 
spin-lattice relaxation times, measured by inversion recovery 
pulse sequence, are (7^,7^) ~ (7.36,10.55) s. The trans¬ 
verse relaxations, obtained by CPMG pulse sequence, have 
characteristic times (7J 1 ,7^) ~ (4.76,0.33) s. We thus study 
processes of maximal duration r ~ 10“ 4 s and consider total 
data acquisition times for the reconstruction of the work and 
entropy statistics within 0.1 ms and 126 ms, being smaller 
than 7^ H,C . This enables us to disregard any energy exchange 
with the system environment during the quenched dynamics. 
The effects of the 13 C transverse relaxation (75p) can be par¬ 
tially overcome by a refocussing strategy in the reconstruction 
procedure. 

We implement the backward process, shown in Fig. 1, by 
driving the system with the time-reversed Hamiltonian, 77® = 
77®_ t , from the equilibrium state, p? — exp(— /377®)/Z r , 
that corresponds to the final Hamiltonian 77® (Z t here de¬ 
notes the partition function at time t). The intermediate 
13 C states during the field quench are p® = U t p^U\ and 
p B _ t = V T -tp Q rV\_ t , where the evolution operators satisfy 
the time-dependent Schrodinger equations d t U t = -m\u t 
and d t Vt = with the initial conditions Uq = Vo = 4. 

Work is performed on the system during the forward and 
backward processes. The corresponding probability distribu¬ 


tions P ¥,B (W) are related via the Tasaki-Crooks fluctuation 
relation [ 29 - 31 ] 

P® (W) /P B (-W) = e ^ w ~ AF) . (3) 

Equation ([3]) characterizes the positive and negative fluctua¬ 
tions of the quantum work W along single realizations. It 
holds for any driving protocol, even beyond the linear re¬ 
sponse regime, and is a generalization of the second law to 
which it reduces on average, (E) = /3((W) — A F) > 0. 

We experimentally verify the arrow of time expressed by 
Eq. 0 by determining both sides of the equation indepen¬ 
dently. We first evaluate the Kullback-Leibler relative entropy 
between forward and backward dynamics by tracking the state 
of the spin-1/2 at any time t with the help of quantum state to¬ 
mography |[23) . Figs.[3](a)-(c) show reconstructed trajectories 
followed by the Bloch vector, for both forward and backward 
processes and different quench times. As a second step, we 
measure the probability distribution P(E) of the irreversible 
entropy production using the Tasaki-Crooks relation Eq. ([3]). 
Employing NMR spectroscopy CD and the method described 
in Refs. |26|-28i (cf. Appendix for a detailed analysis), we 
determine the forward and backward work distributions, from 
which we extract (3, W and AP, and hence the entropy pro¬ 
duced during each process. The measured nonequilibrium en¬ 
tropy distribution is shown in Fig. [4] (a). It is discrete as ex¬ 
pected for a quantum system. We further observe that both 
positive and negative values occur owing to the stochastic na¬ 
ture of the problem. However, the mean entropy production 
is positive (red line) in agreement with Clausius inequality 
(E) > 0 for an isolated system. We have thus directly tested 
one of the fundamental expressions of the second law of ther¬ 
modynamics at the level of an isolated quantum system (3) . 

A comparison of the mean entropy production with the 
Kullback-Leibler relative entropy between forward and back- 
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ward states is displayed in Fig. 0(b) as a function of the 
quench time. We observe good agreement between the two 
quantities within experimental errors that are due to inhomo¬ 
geneities in the driving rf-field and non-idealities of the field 
modulation. These results provide a first experimental con¬ 
firmation of Eq. (1) and the direct verification of the thermo¬ 
dynamic arrow of time in a driven quantum system. They 
quantify in a precise manner the intuitive notion that the faster 
a system is driven away from thermal equilibrium (i.e. the 
bigger the mean entropy production or the shorter the driving 
time r), the larger the degree of irreversibility, as measured by 
the relative entropy between a process and its time reversal. 

In the linear response regime m, Onsager has derived 
generic expressions for the entropy production which form 
the backbone of standard non equilibrium thermodynamics. 
These results are, however, limited to systems that are driven 
close to thermal equilibrium. By contrast, Eq. 0 holds for 
any driving protocol and thus arbitrarily far from equilibrium. 
In order to check the general validity of Eq. 0, we use the 
linear response (LR) approximation of the mean work na, 
(Wul) = A F + (3AW 2 /2, where A W 2 is the variance of 
the work, to obtain the mean entropy production (Elr) = 
/3 2 AW 2 /2. Fig. 0(c) shows the experimental values of (E) 
and (E L r) as a function of the quench duration. We note 
that the measured irreversible entropy production (E) is close 
yet systematically distinct from its linear response approxi¬ 
mation (Elr), the difference being more pronounced for fast 
quenches (small r), as expected. Fig. [4] (c) thus suggests that 
the quenches implemented in the experiment are performed 
somewhat away from the nonlinear response regime. We also 
mention that we achieve good agreement between experimen¬ 
tal data (dots) and numerical simulations (dashed lines). 
Conclusions.- We have assessed the emergence of the arrow 
of time in a thermodynamically irreversible process by using 
the tools provided by the framework of non equilibrium quan¬ 
tum thermodynamics. We have implemented a fast quenched 
dynamics on an effective qubit in an NMR setting, assessing 
both the mean entropy produced across the process and the 
distance between the state of the system and its reverse ver¬ 
sion, at all times of the evolution. Let us discuss the physical 
origin of such time-asymmetry in a closed quantum system. 
Using an argument put forward by Loschmidt in the clas¬ 
sical context, its time evolution should in principle be fully 
reversible m. How can then a unitary equation, like the 
Schrodinger equation, lead to Eq. 0 that contains a strictly 
nonnegative relative entropy? The answer to this puzzling 
question lies in the observation that the description of phys¬ 
ical processes requires both equations of motion and initial 
conditions mini. The choice of an initial thermal equilib¬ 
rium state singles out a particular value of the entropy, breaks 
time-reversal invariance and thus leads to the arrow of time. 
The dynamics can only be fully reversible for a genuine equi¬ 
librium process for which the entropy production vanishes at 
all times. Moreover, issues linked to the "complexity" of the 
preparation of the initial state to be used in the forward dy¬ 
namics (or the corresponding one associated with the time- 


reversed evolution) have to be considered f32l . By provid¬ 
ing an experimental assessment of the microscopic foundation 
of irreversibility (systematically beyond the linear response 
regime), our results both elucidate and quantify the physical 
origin of the arrow of time in the quantum setting of an iso¬ 
lated system. 

APPENDIX 

1. Initial thermal state preparation 

Preparation of the initial state can be done employing spa¬ 
tial averaging techniques, i.e., a combination of transverse 
rf pulses and longitudinal field gradients [23, 33 ]. Through 
this standard method, we prepared a joint pseudo pure state 
equivalent to |0) H (0| (8) for the Hydrogen and the Carbon 
nuclear spins, with a = 0, r for the forward and backward 
protocols, respectively. We are using the logical representa¬ 
tion of the ground (|0)) and excited (|1)) spin states. The 
Hydrogen nuclear spin is prepared in the ground state while 
the Carbon is effectively prepared in a thermal diagonal state 
p°a = e~P n <x / Z a , where Z a = tre - ^" is the partition 
function for . The initial populations of in the ba¬ 
sis follow the Gibbs distribution associated with inverse spin 
(pseudo-)temperature /3 = (&bT) _1 . 

In all experiments we fixed the spin temperature of the 
Carbon in the initial Gibbs thermal distribution, such that 
ksT/h = 1.56 ± 0.07 kHz. The uncertainty in this value 
is obtained through a comparison of Carbon populations in 
a series of independent state preparation and quantum state 
tomography. The trace distance between the experimentally 
prepared initial state (p exp ) and the ideal Gibbs state (p ldeal ), 
defined as S = tr |p exp — p ldeal |i/2, is smaller than 0.03 for 
all prepared states, in other words a less than 3% chance to 
discriminate them. 

2. Work and irreversible entropy distributions. 

The work performed on or by a closed quantum system un¬ 
dergoing a protocol (dictated by ldF t ^) is a stochastic vari¬ 
able ifTOl [Till . For a quasi-static and isothermal process, the 
average work equals the variation in free energy. On the other 
hand, in a scenario involving a non-quasistatic and driven 
closed system dynamics (as in our experiment), an amount 
of work may be done on the system above the variation in 
free energy. This introduces a dissipated work , (Wdi SS ) = 
(W) — A F > 0, inherently linked to the irreversible nature 
of the process. Such irreversibility is associated to an entropy 
production 18UT3L 

E = p (W - A F) . (4) 

which is also a stochastic variable that depends on the work 
distribution, the net change in the system free energy, and the 
inverse temperature. So the entropy production probability 
distribution, P(E), is built upon the experimental assessment 
of all the mentioned variables. 
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Figure 5. NMR pulse-sequence for the reconstruction of the work characteristic functions x F (^) and x B (u)- a, Sequence for the forward 
process (a — 0). b, Sequence for the backward process (a — r). We start from the pseudo-pure state |0) H (0| (8) p *£. The blue (red) circles 
represent transverse rf-pulses in the x ( y ) direction that produce rotations by the displayed angle. The orange connections represent free 
evolutions under the scalar interaction, l~ij — 2nJcr^a^ (with J « 215.1 Hz), during the time displayed above each connection. This 
time-length of the coupling is varied as function of s, which is related to the conjugate variable u in Eq. (5]) as s — 2 'kv$ u. 


An important step in our experiment is the characterization 
of the probability distribution of work fluctuations during the 
quenched dynamics. Such a distribution is obtained using a 
Ramsey-like interferometric method theoretically proposed in 
refs. ED 128) and experimentally implemented to verify the 
fluctuation relations in a quantum scenario f26l . This ap¬ 
proach is based on the interferometric measurement of the 
characteristic function of the work distribution 134) . which 
can be written for the forward process (F) as 

= ( 5 ) 

m,n 

where pn q, ° = e _/3en /Z 0 is the initial occupation probabil¬ 
ity for the n-th eigenstate (|no)) of H ¥ with energy e n , while 
P F m\n = l( m r|^r|^o) | 2 is the conditional transition proba¬ 
bility to find the system in the ra-th eigenstate (| m T )) of H ¥ 
(with energy e m ) if it was previously in |no) at t = 0. 

The characteristic function can also be written as 

X F (^) = tr [(U T e- iun o)p^(e- iu <U T )^ . (6) 

The main idea of the protocol to reconstruct x F (u) is to es¬ 
timate the trace in Eq. via an ancillary system ll27l 28l. 
the Hydrogen nuclear spin in the present experiment. The 
NMR pulse sequence employed to perform such a task is de¬ 
picted in Fig. [5] The ancillary system, initially prepared in the 
ground pseudo-pure state, is rotated to an equal-weighted su¬ 
perposition (equivalent to (|0) H + |1) H ) /y/2) in the grey box 
of Fig. [5] The Carbon nuclear spin is the driven system of 
our interest. For technical reasons (see Sect. 1 above), the 
13 C spin starts from a mixture of computational-basis (a^- 
basis) states, with occupation probabilities given by pn q ’° (or 
by Pn q,r in backward mode). In fact, the state basis will be fur¬ 
ther rotated to the initial Hamiltonian one. Controlled-unitary 


operations are performed using the natural scalar coupling in¬ 
teraction % j = 2tt Jg^g^, which is proportional to g^ , while 
the initial and final Hamiltonians H ¥ and H ¥ in the forward 
protocol are proportional to Gy and cr£, respectively. This is 
compensated by the rotations inside the green boxes of Fig. [5] 
Such rotations also account for the fact that the initial state 
was prepared in the -basis. The same reasoning is applied 
in the backward case. The purple boxes in Fig. [5] are a refo¬ 
cus strategy to mitigate the effects of the transverse relaxation. 
The quenched dynamics is a consequence of a suitable mod¬ 
ulation (amplitude and phase) of a transverse rf field, which 
can be described by the Hamiltonian in Eq. (1) for the for¬ 
ward protocol. The final step of the algorithm is the measure¬ 
ment of the free induced decay (FID) signal of the Hydrogen 
nuclear spin. From this signal one can obtain the transverse 
magnetization, where the characteristic function is encoded 
as ^f,b ^ _ 2 2 1 (o^y Application of an inverse 

Fourier transform allows us to obtain the work distribution for 
the quenched dynamics. 

We have measured several experimental configurations, 
keeping fixed the initial spin temperature (given by the 
weights Pn q,0( ^ of the initial Carbon population), and varying 
the quench type (forward or backward) and the quench du¬ 
ration. For each configuration, the interaction time s of the 
free evolution under the scalar coupling in Fig. [5] was var¬ 
ied through 360 equally-spaced values; each realisation cor¬ 
responds to an independent experiment with an average over 
a set of identical and independent molecules. 

A typical output of the characteristic function reconstruc¬ 
tion algorithm is shown in Fig. [6] (a), where each data point 
corresponds to an independent experiment for each scalar 
interaction parameter s. The inverse Fourier transform of 
the transverse magnetization (the work distribution P F (W)) 
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Figure 6. Typical output of the interferometric circuit shown in 
Fig. [5] for a forward quench of time-length r = 100 /is. a, Blue 
and red symbols (lines) correspond to the real and imaginary parts of 
the characteristic function x F (u) (where u = (27ri/o) _1 s), i.e., the 
experimentally measured x and y components of the 1 H transverse 
magnetization displayed in arbitrary units (a. u.), as function of the 
adimensional parameter s. b, Work distribution P F (W) obtained 
from the inverse Fourier transform of the the characteristic function, 
c and d, Same as panels a and b, but after subtracting the noise in the 
Fourier analyses. 

is displayed in Fig. 0(b). From Eq. ([5]), we would ex¬ 
pect to observe four peaks, since we have a quenched two- 
level system. Furthermore, we expect the peaks to be lo¬ 
cated around {±(z/o ± z+)}> which correspond to the fre¬ 
quencies {—2.8, —0.8, +0.8, +2.8} kHz considering the pa¬ 
rameters used in our setup (zx 0 = 1.0 kHz and u T = 1.8 kHz). 
The peaks in Fig. [6] (b) are proportional to the probabilities 

Pi q ’°Po\y Po q ’Vojo> Pi+ijk and Po+ijo’ from left t0 ri g ht > 

respectively. 

A least-squares method can be applied to estimate the 
amplitudes of the four peaks in Fig. [6] (b), which can be 
identified with the terms Pn q ’ 0 p^q n of Eq. (pi). We can also use 
this data to filter the noise frequencies in the original char¬ 
acteristic function (cf. Fig. [6] (a)] disregarding contributions 
of other frequencies. The noise sources are non-idealities in 
the refocus sequence of the scalar interaction (purple boxes 
in Fig. 0). The noise-subtracted characteristic function and 
work probability distribution are shown in Figs. |6j: and[6ji, 
respectively. The typical noise amplitude was added to the 
error bar in the estimation of the total probabilities Pn q ’°P^[ n - 

3. Estimation of the net difference in free energy. 

The Crooks relation can be used to experimentally estimate 
the free energy difference in our nonequilibrium dynam¬ 
ics. We have performed a set of experiments for different 
quenches with durations of r = 100, 200, 260, 320, 420, 500, 
and TOO /xs. For each quench duration, we performed both the 


Figure 7. Estimation of the net difference in free energy. Crooks rela¬ 
tion P F (+f+)/P B (— W) (in logarithmic scale) for different quench 
time-lengths. Symbols represent experimental data; the dashed line 
and shaded regions represent the best linear fit and its estimated un¬ 
certainty, respectively. The net difference in the free energy, AP, 
can be obtained from the interception of the curves and the yellow 
dashed line (where P F (W) m P B (-W) and (3W = j3AF). 

forward and backward protocols (as described in Fig. [5]). The 
ratio between the forward and backward work distributions 
are plotted in a logarithmic scale in Fig. [7] for the different 
quenches, where we found the best fitting line to the equation 
(which represents the Crooks relation): 

"• j+r= m - af > ■ <7) 

The data points in Fig. [7] represent the ratio 
P F (+VE)/ P B (-W) only for the inner peaks of the work dis¬ 
tributions. Since these points have themselves an uncertainty 
(from the Fourier analyses), the linear adjustment also has an 
uncertainty, which is represented by shaded regions in Fig. [7] 
The angular coefficient of the linear fitting is the inverse spin 
temperature [3 while the linear coefficient is the product (3 AF. 

4. Entropy production expanded in terms of the cumu- 
lants of the work distribution. 

The Jarzynski equality allows one to write the free-energy dif¬ 
ference as AF = —/7 _1 ln(e _ ^ w/ ), which can be plugged in 
Eq. ([4]) resulting in the following expansion 

(£) = E M/5)/?" ’ (8) 

n> 2 

with { ft n (/3)} the set of cumulants of the work distribu¬ 
tion. The first two cumulants are the work average, ft 1 = 
(W), and the work-distribution variance, ft 2 = (AW) 2 = 

(W — (W)) 2 \ respectively. As a first approximation, 
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the mean produced entropy is proportional to the work- 
distribution variance. If the system presents a linear response 
to the driving protocol, the relation (Elr) = P 2 (A W) 2 /2 
turns out to be a very good approximation. Moreover, at 
macroscopic scales, the total work involves, in general, a sum 
of many independent random variables. In this case, the Cen¬ 
tral Limit Theorem implies that the total work distribution 
converges to a Gaussian one, for which all cumulants of or¬ 
der higher than two vanish, making the above approximation 
exact. 

In our experiment, we are dealing with the microscopic be¬ 
haviour of a single quantum system driven by a quenched evo¬ 
lution, which leads to the entropy production distribution (and 
also the work distribution) notably non-Gaussian, as can be 
observed in Fig. 4 (a) of the main text and also in Fig. [6] Fur¬ 
thermore our spin system responds in a non-linear way to the 
field quench as noted in Fig. 3 (a), (b) of the main text. The 
system magnetization acquires a longitudinal component dur¬ 
ing the dynamics driven by a fast time-varying Hamiltonian 
in the transverse direction. These facts introduce a finite de¬ 
viation from the expectations of linear response theory as ob¬ 
served in Fig. 4 (c) of the main text. We can also note in 
this figure that for a “long duration” quench the mean entropy 
produced tend to a finite value, since we have a closed system 
dynamics (it has no thermalisation along the quench protocol). 

The numerical simulations curves, dashed lines plotted 
in Fig. 4 (c) of the main text, were obtained through a nu¬ 
merical integration of the conditional transition probabilities, 
P F rn\n = l( m r|^r|^o) | 2 (employing a code written in Python 
with the aid of the QuTip toolbox [35]), from which we can 
compute P(W ), P(E), (E), and (T, LR ). 

5. Derivation of the microscopic expression for the entropy 
production (Eq. (1) of the main text). 

For the sake of completeness we will discuss here a derivation 
of Eq. (1) in the main text, which will also be very instructive 
for the purpose of the present article. Let us consider the sce¬ 
nario described in Fig. 2 of the main text. The mean entropy 
production (E) = ft (W) — /3AF in the forward protocol can 
be written as 


This derivation did not assume any particular relationship be¬ 
tween p eq and p ¥ r . The assumption of a closed, unitary dy¬ 
namics implies that von Neumann entropies of initial and fi¬ 
nal states are equal, and the irreversible entropy produced dur¬ 
ing the quench is exactly given by the Kullback-Leibler diver¬ 
gence of the final nonequilibrium state and the equilibrium 
state of the final Hamiltonian. 

In addition, the assumption of unitary dynamics allows us 
to prove that the Kullback-Leibler divergence between the for¬ 
ward and backward dynamics is constant during the quench, 
i.e., 5 (^|| pr q ) = S (p?|| with 

S (p F || p*_ t ) = tr (p F In p F ) - tr (p F In p*_ t ) . (12) 

We note that the intermediate states during the quench are 
p v t = U t p q U\ and p®_ t = V T - t p e Pvl_ t (for the forward and 
backward protocol, respectively), where the evolution oper¬ 
ators satisfy d t U t = —ibi-Mt and d t Vt = iH, r -tVt (with the 
initial conditions Uo = Vo = H). A simple change of variables 
leads to d t V T - t = — il~i t V T -u which is the same equation 
as that for U t , but with a different initial condition, namely 
that V T - t | t=r= H- By inspection of the structure of the 
Schrodinger equation, we are allowed to right-multiply one 
solution by any constant operator to obtain another solution 
for a different initial condition. Doing that, we obtain a rela¬ 
tion between the forward and backward dynamics as 

V T -t = U t U\ . (13) 


To prove Eq. (1) of the main text, we shall investigate the time 
evolution of both terms in Eq. (12). The term tr (/ p ¥ t In p f) is 
constant because of the unitary evolution. It remains to show 
that the term tr (p* In p^_ t ) is also constant, which can be 
proven as follows 


tr [pt InPr-t] 


tr 

tr 


Utp^Ul In (v T - t p e M_ t ) 
fffl4v T -t In (p e T q ) V\_ t U t 


tr [p^U\ In (/ o f?) Ur] 

tr [Pr In (Pr q )] , (14) 


(S> = p tr (n F T p F ) - (3 tr (?«) - /3AF , (9) 

where AF = /3 -1 In (Zo/Z T ) is the free energy differ¬ 
ence between the equilibrium states p sq = e _/3 ^°/Zo and 
Pr q = e~ f3n ^ /Z r , at constant inverse temperature p. We can 
identify 

In p e t q — U/3 -1 In Z t , (10) 

for t = 0 and t = r. Using the identification ( fT0| ) in Eq. 0. 
we obtain (E) = — tr (p^ In p^ q ) + tr (pg q In pj^J, which can 
be written as a sum of the Kullback-Leibler divergence and a 
variation in von Neumann entropy [<S V n (p) = — tr(p In p)], 

<s> = <5 (p F T \\p e T q ) + <S vN (^) - <S vN (p e o q ) . (11) 


where Eq. ( [13] ) was used in the last-but-one equality. This 
shows that <Sl[p£|| p®_ t ) is time invariant. In this way we have 
verified Eq. (1) in the main text. 

For a prior derivation of this result, employing different ar¬ 
guments and methods, see Ref. iTbl . 
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